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Abstract 

We present a new global optical potential (GOP) for nucleus-nucleus systems, including neutron-rich and 
proton-rich isotopes, in the energy range of 50 ~ 400 MeV/u. The GOP is derived from the microscopic 
folding model with the complex G-matrix interaction CEG07 and the global density presented by Sao Paulo 
group. The folding model well accounts for realistic complex optical potentials of nucleus-nucleus systems 
and reproduces the existing elastic scattering data for stable heavy-ion projectiles at incident energies above 
50 MeV/u. We then calculate the folding-model potentials (FMPs) for projectiles of even-even isotopes, 
^-22c, 12-24q^ i6-38Ne, 20-40]y[g^ 22-483j^ 26-525^ 30-62^j.^ 34-70(~.^^ Scattered by stable target nuclei 
of 12c, 1^0, 28si, 40ca 58Ni, 90zr, 1203^^ ^nd ^o^pb at the incident energy of 50, 60, 70, 80, 100, 120, 
140, 160, 180, 200, 250, 300, 350, and 400 MeV/u. The calculated FMP is represented, with a sufficient 
accuracy, by a linear combination of 10-range Gaussian functions. The expansion coefficients depend on the 
incident energy, the projectile and target mass numbers and the projectile atomic number, while the range 
parameters are taken to depend only on the projectile and target mass numbers. The adequate mass region 
of the present GOP by the global density is inspected in comparison with FMP by realistic density. The full 
set of the range parameters and the coefficients for all the projectile-target combinations at each incident 
energy are provided on a permanent open-access website together with a Fortran program for calculating 
the microscopic-basis GOP (MGOP) for a desired projectile nucleus by the spline interpolation over the 
incident energy and the target mass number. 

PACS numbers: 24.50.-Hg, 24.10.Ht, 25.70.Bc, 25.70.-z 
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I. INTRODUCTION 



The elastic scattering process contains valuable information about nuclear many-body dynam- 
ics in nuclear reactions induced by a nucleon or a composite-nucleus projectile on a target nucleus. 
The elastic scattering observables (averaged over a certain range of incident energy) are known to 
be well described by a local or non-local one-body type complex potential called the optical po- 
tential. The shape and the strength of the complex optical potential and their dependence on the 
incident energy and the mass and charge of the colliding system reflect, on one hand, the static 
properties of nuclear structure of the colliding system and the basic nucleon-nucleon (NN) inter- 
action and, on the other hand, the effects of complicated non-elastic reaction processes that follow 
the elastic scattering process [|i|-|4|]. 

The optical potential is also indispensable to a reliable extraction of nuclear structure infor- 
mation from experimental data of various direct nuclear reactions such as inelastic scattering 
and transfer reactions. Namely, most direct reaction theories, such as the distorted-wave Bom 
approximation (DWBA), distorted-wave impulse approximation (DWIA) and coupled-channels 
(CC) method, incorporate the optical potential as the distorting potentials in the entrance or exit 
channels in the analyses of experimental data M\- It is often the case that the calculated results 
strongly depend on the choice of the distorting potential. Therefore, it is very important to es- 
tablish reliable optical potentials (or more generally interaction models between colliding nuclear 
systems) to survey unknown properties of nuclear structure and reaction dynamics through the 
experimental data. 

Historically, phenomenological optical potentials (POPs) were deduced by analyzing individual 
elastic scattering data for various targets and incident energies. The systematic analyses of a large 
body of scattering data enabled us to find systematic behavior of potential parameters with the 
variation of the mass and charge of the target nucleus as well as the incident energy, which led 
to a proposal of the global optical potential (GOP). The GOP for nucleon-nucleus scattering has 
a long history for more than four decades [|5hZ|] and nowadays it is in almost established stage 
covering whole stable target nuclei and incident energies from nearly zero to 1 GeV [|8l-ll2||. For 
deuteron-nucleus scattering, various GOPs have been proposed [13-[l7] that are also applicable to 
a wide range of incident energy and target mass number. 

The study of optical potentials for composite projectiles heavier than deuteron is still in a devel- 
oping stage and far from satisfactory level, partly because of insuflicient number of experimental 
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data on elastic scattering compared with nucleon and deuteron cases and, more essentially, be- 
cause of the composite nature of the system which makes the nuclear reaction dynamics more 
complicated, such as mutual excitations, breakup, transfer, fusion etc. The increase of the number 
of reaction channels leads to a strong absorption of flux from the elastic-scattering channel, which 
makes it diflicult to uniquely determine the shape and strength of the optical potential yD. 

Nevertheless, a number of eff"orts have been made to deduce the optical potentials for com- 
posite projectiles heavier than deuteron and GOPs were also proposed for some light ions, such 
as ^He OH, a particles [|20l- l22ll and some light heavy-ion projectiles ^'^Li [|23n and ^He 11241 1. 
However, the range of the incident energy and target mass number covered by those GOPs is very 
limited compared with those for nucleon and deuteron scattering. 

For heavy-ion projectiles, a number of elastic- scattering experiments have been done and op- 
tical potentials have been deduced where ambiguity and uniqueness of the deduced optical po- 
tentials were discussed and some sort of systematics were studied n25\\ . However, the number of 
experimental data is far from satisfactory level to extract GOPs, except for some limited scattering 



systems for which the proposed "GOPs" are only applicable to a limited energy range [|26 



Nowadays the experimental facilities have been developed intensively. The radioactive isotope 
beam is produced with a considerable intensity up to the driplines for light elements and its energy 
reaches up to 400 MeV/u. In order to deduce meaningful information on the properties and nuclear 
structure of such rare isotopes through reaction experiments, it is essentially important to establish 
the optical potentials, or any alternative interaction models, for those unstable heavy-ion systems 
up to the highest energy region, say 400 MeV/u. In most radioactive ion beam experiments, 
however, it is very difficult to measure reliable elastic-scattering cross sections, particularly to 
obtain a precise angular distribution necessary for extracting reliable optical potentials. Therefore, 
one needs a reliable GOP or any alternative theoretical interaction models. 

The global description of the optical potential has been attempted in the view point of the 
microscopic interaction models. For a nucleon-nucleus system, the most successful description 
of the complex optical potentials and the elastic scattering data is obtained by the single-folding 
model approach with complex G-matrix interactions [|28l-l39|l based on the Bruckner Hartree-Fock 
theory. The energy and target-mass dependences of the geometrical shape and strength of the 
potential are explained by various kinds of medium eff"ects of effective interaction (G-matrix) in 
nuclear medium as well as the exchange terms and the non-locality, in addition to the genuine 
energy dependence of the free-space A'^A'^ interaction used and the nucleon density distribution of 
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the target nucleus. 

Microscopic description of the nucleus-nucleus interaction has also been studied for decades 
on the basis of double-folding model (DFM) with the G-matrix interactions. The most successful 
models in the earlier stage are those based on the M3Y interaction [4^, 41] or its extension having 



various types of density dependence nlSi 
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42I-I45II. However, all these interaction models for 



heavy-ion optical potential predict only the real part of the potential and the phenomenological 



imaginary part has to be added by hand. The JLM interaction [13 311 had been the only comp 



matrix interaction applied to the construction of complex heavy-ion optical potential ll46l - l49ll . The 



ex G- 



JLM interaction is known to be successful for nucleon-nucleus system but it has a crucial problem 
in the application to the nucleus-nucleus systems. Namely, the calculated DFM potential with the 
JLM interaction is too strong both in the real and imaginary parts compared to the POP evaluated 
from experimental data and needs a strong renormalization. The origin of the renormalization 
problem in the DFM with JLM interaction is found to be originated from impro per p rescriptions 
(average prescriptions) for evaluating the local density in the DFM calculation [49], which was 
inevitable for the JLM interaction that was given only for nuclear matter with density below the 
saturation density po. 

In order to overcome this problem, three of thepresent authors have recently constructed a 



new complex G-matrix interaction called CEG07 [|38l Is^, that is applicable for nucleon density 
up to twice the saturation density, p ~ 2po in the energy range of £ = 20 ~ 400 MeV/u and 
they have proposed the DFM with the CEG07 interaction that gives a reliable complex optical 
potential for nucleus -nucleus systems [S^S^- It is shown that the DFM with the new G-matrix 
predicts the proper strength and the shape of the nucleus-nucleus optical potentials and reproduces 
almost all the existing experimental data for elastic scattering of '^C and '^0 projectiles by var- 
ious stable target nuclei in the energy-range of £ = 70 ~ 135 MeV/u HE 1|] with essentially 
no renormalization. It is interesting to note that the folding model also predicts that the real part 
of heavy-ion optical potentials changes its sign from negative (attractive) to positive (repulsive) 
around E = 200 ~ 300 MeV/u [|52|l . In addition to the elastic scattering, the CEG07 folding model 
has also been applied to the coupled-channel analyses of inelastic scattering with considerable suc- 



cess [1530. The present DFM with CEG07 is now the most reliable microscopic interaction model 
at present for constructing the complex optical potential for nucleus-nucleus scattering systems. 

However, one needs rather complicated and elaborated numerical tasks to calculate the DFM 
potential with CEG07 because of its precise treatment of density dependence as well as that of the 
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FIG. 1 : (Color online) The table on the web site 115411 . 



finite-range exchange component and, hence, it will not be a better choice to perform heavy numer- 
ical calculations to analyze individual scattering experiments for each combination of projectile- 
target system at each scattering energy. It would rather be desirable to construct reliable GOP on 
the basis of the reliable microscopic interaction model. 

In the present paper, we apply the present microscopic folding model with the CEG07 interac- 
tion to the construction of GOP for scattering of heavy ions including neutron-rich and proton-rich 
unstable isotopes over the wide range of nuclear chart. To this end, we first calculate the DFM 
potentials for projectiles of even-even isotopes, ^"^^C, ^^'^"^O, ^^"^^Ne, ^"""^"Mg, ^^"'^^Si, ^''"^^S, 
3()-62^j.^ and 34-70(-a (see Fig. [B, scattered by stable target nuclei of ^^C, ^^O, ^^Si, ^"Ca ^^Ni, 
^^Zr, i^t'Sn, and ^o^Pb at certain steps of incident energy, say E = 50, 60, 70, 80, 100, 120, 140, 
160, 180, 200, 250, 300, 350, and 400 MeV/u. The calculated folding potential is represented, 
with a sufficient accuracy, by a linear combination of 10-range Gaussian functions. The expan- 
sion coefficients depend on the incident energy, the projectile and target mass numbers and the 
projectile atomic number, while the Gauss-range parameters are taken to depend only on the pro- 
jectile and target mass numbers. The full set of the range parameters and the coefficients for all the 
projectile-target combinations at each incident energy are then stored on a permanent open-access 
website [|54|] together with a Fortran program for calculating the microscopic -basis GOP (MGOP) 
for a requested projectile nucleus by the spline interpolation of the parameter sets over the incident 
energy and the target mass number. These parameter sets and the program will be provided to open 
use for a wide range of nuclear reaction analyses. 

In the next section, we briefly describe the theoretical frame of the microscopic folding model 
with the CEG07 interaction together with the parameterizations of the calculated DFM poten- 
tials in a form of GOP, where we discuss the global parametrization of the input nucleon density 
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distributions. In sect.III, we first compare elastic scattering cross sections calculated with the mi- 
croscopic DFM potential as well as its equivalent MGOP with the existing experimental data for 
stable-nucleus heavy-ion projectiles, ^^C, ^^O and "^"Ar, by various target nuclei at various incident 
energy, to confirm the applicability of the present MGOPs. We them calculate elastic- scattering 
cross sections with the MGOPs for projectile of unstable nuclei including those near the driplines 
and discuss the isotope dependence of the predicted cross sections that reflect difi'erent density 
distributions through the DFM procedure. The final section will be devoted to conclusion. 



II. GLOBAL OPTICAL POTENTIAL 



A. Double folding model 

We construct the nucleus-nucleus optical model potential based on the DFM using the complex 



G-matrix interaction CEG07. In the previous work [|50Ld1|| . it was found that the three-body force 
effect, particularly its repulsive component, plays a crucial role to obtain the proper strength and 
shape of the nucleus-nucleus potential consistent with the experimental data. Thus, we use the 
CEGOVb version of the CEG07 interaction that includes the three-body force effect. 

The microscopic nucleus-nucleus potential can be written as a Hartree-Fock type potential: 

Uv = X '^^ OIvdIO" > + < OIvexI;? >] (1) 

ieAi,jeA2 

= f/D + C/ex, (2) 

where Vd and vex are the direct and exchange parts of complex G-matrix interaction. The exchange 
part is a nonlocal potential in general but can be treated as a local potential by using the plane- wave 



representation for the A^A^ relative motion 115 5l 15611 . The direct and exchange parts of the localized 



potential are then written in the standard form of the DFM potential as 

Uv,{R) = J pi{ri)p2{r2)vu(s;p,E)dridr2 (3) 

= J [pf\r,)pf\r2Wf(s;p,E) 
+pf{r,)pf{r2)vf\s;p,E) 
^pf{r,)pf{r2)vt'\s;p,E) 

+d^\r,)pf{r2)v'^'\s-p, E)\dr,dr2, (4) 
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where s = r2 - Vi + R, and 
UExiR) 



J Pi(ri,ri + s)p2ir2,r2 - s)vEx(y,P,E) 
ik(R) ■ s 



X exp 



M 



dridr2 



(5) 



pf^Cn, n + s)p'^>{r2, r2 - sy^^\s;p, E) 
+pf\rur, + s)pf(r2, r2 - s)vS\^;p, E) 
V/'^Cr-i, ri + s)pf(r2, r2 - s)v^^^\s-p, E) 
+pf\ri,ri + s)pf \r2, r2 - s)v^^^\s;p, £)} 



xexp 



ik(R) ■ s 
M 



dridr2. 



(6) 



where E is the incident energy per nucleon (MeV/u). The superscript of p (p and n) indicates 
the proton and neutron densities, respectively. The direct (D) and exchange (EX) parts of the 
proton-proton (pp), proton-neutron (pn), neutron-proton (np), and neutron-neutron (nn) G-matrix 
interactions are written as 



(pp,nn) _ 



1 



D,EX 



(7) 
(8) 



= 4(v°'±3v'^), 

= ^(±v°° + v°>+3v'°±3v"), 
8 

in terms of v^^ , the spin-isospin component (5=0 or 1 and T=0 or 1) of the G-matrix interaction. 
Here, the upper and lower part of the double-sign symbols (+) in the r.h.s. correspond to the 
direct (D) and exchange (EX) parts, respectively. Here, the exponential function in Eq. ^ is 



transformed into the spherical Bessel function, 7o(^^), by the multipole expansion. The k(K) is 



the local momentum for nucleus -nucleus relative motion, the magnitude of which is given by 



k\R) = 



2mM. 



-[£c.m.-Ret/F(i?)-Koui(^)], 



(9) 



where M = AiA2/(Ai -i- A2), ^c.m. is the center-of-mass energy, m is the nucleon mass and Vcoui 
is the Coulomb potential. Here, the Coulomb potential Vcoui is also obtained by folding the A^A^ 
Coulomb potential with the proton density distributions of the projectile and target nuclei. Ai and 
A2 are the mass numbers of the projectile and target, respectively. The exchange part is calculated 
self-consistently on the basis of the local energy approximation through Eq. The density 
matrix p(r, r') is approximated in the same manner as in Ref. nSln : 



P(-^^'^ = ^Mkfs)p(^^), 



(10) 



where kSj^ is the effective Fermi momentum [15 811 defined by 



,eff 



3p2 



(11) 



where we adopt Cs = 1/4 following Ref. |!590. In the present calculations, we employ the so-called 
frozen density approximation (FDA) for evaluating the local density. In the FDA, the density- 
dependent NN interaction is assumed to feel the local density defined as the sum of densities in 
colliding nuclei evaluated at the mid-point of the interacting nucleon pair: 



P= Pi+Pi 



(12) 



where the local densities are evaluated at the position of each nucleon for the direct part and at 



the middle point of the interacting nucleon pair for the 



been widely used also in the standard DFM calculations [|4lL |44L 



exch ange part [160 ^ 



45 



50 



59 



m. The FDA has 



611]. In Ref. [50J it is 



confirmed that FDA is the best prescription in the case with CEGOVb to reproduce the data. 



B. Input density 

The input density is important to construct the global optical potential in the folding procedure 



In this paper, we adopt the Sao Paulo density provided by Sao Paulo group [|26ll as a global density. 



The Sao Paulo density has the functional form of the two-parameter Fermi, as 

Pv,n{r) = ^° > (13) 

l+exp(^) 

where, the subscript of p, R, and a (p and n) indicates the proton and neutron parts, respectively. 



The parameters of i?p,n and obtained in Ref. [|26|] . as 



7?p = 1.81Z^^^ - 1.12, (14) 

i?„ = 1.49iV^^^ -0.79, (15) 

ap = 0.47 - 0.00083Z, (16) 

an = 0.47 + 0.00046A^, (17) 

where, Z and A'^ are the proton and neutron number, respectively. 

Here, we note that the Sao Paulo density is fitted to the calculated results of the Dirac-Hartree- 
Bogoliubov model in the stable nuclear region. It is nontrivial that the Sao Paulo density is useful 



in the unstable nuclear region. However, our purpose is the construction of global optical potential 
with various projectiles including the proton-rich or neutron-rich nucleus. Therefore, we compare 
the results obtained by the Sao Paulo density with those by a theoretical framework aimed at 
precise description of nucleon density for unstable nucleus region up to driplines. As shown in 
Sec. III. B, the GOP with the Sao Paulo density gives the similar results to the results by the FMP 
with such realistic density, except in the extreme case, such as near or on the dripline. For such 



64|], for 



realistic densities, we adopt those proposed by Niigata group for Carbon isotopes 1162 - 
Oxygen isotopes n65\\ . and for ^"Ne nucleus ll66ll . and we call them as Niigata density. Niigata 
densities have already been applied to the analyses of total reaction cross sections in the Glauber 
model and found to well reproduce the experimental data [1641 16511 . The Niigata density is obtained 
by a Slater determinant based on empirical nucleon separation energies. A core -i- two-nucleon 
model is employed for ^^C [62] and ^^C [63I in order to take into account some dynamical effects 
beyond the simple mean field model. Since ^^C is unstable with respect to a neutron emission, it 
is more realistic to describe ^^C with the three-body model. 

Figure |2] shows the proton, neutron, and matter density distributions of the Sao Paulo and 
Niigata densities for the Oxygen isotopes. The two kinds of density distributions have similar form 
near the stable region, especially for the matter distribution. On the other hand, the difference is 
clearly seen on the dripline. The large difference in the middle and long range parts of the matter 
density distribution reflects the difference in the neutron density distribution. In the next section, 
we compare the Sao Paulo density with the Niigata one in the calculated elastic cross section and 
investigate the adequate mass region of the present GOP with the Sao Paulo density. 



C. Functional form of global optical potential 

We parameterize the nuclear part of the DFM potential in a functional form. The real and 
imaginary parts of the calculated DFM potential are represented in terms of a linear combination 
of 10-range Gaussian form, respectively, with sufficient accuracy as 

exp(-^)} =Vgop(R), (18) 



( 

n=l ^ 

WAR) = Yj {fi" (-^)} = ^GOp(^), (19) 



10 




0246 80 246 80 2468 
r(fm) 



FIG. 2: Comparison the Sao Paulo (SP) density with the Niigata one for the Oxygen isotopes. The soUd and 
dotted curves are the Niigata and Sao Paulo (SP) densities, respectively. The alphabets of the subscription 
(a, b, and c) indicate the proton, neutron, and matter densities, respectively. The numbers of the subscription 
(1, 2, 3, and 4) indicate the mass number {A = 18, 20, 22, and 24) for the Oxygen isotopes, respectively. 



where, Vf(^) and Wp(R) are the real and imaginary parts of the FMP. The Vgop(R) and Wqop(R) 
are defined as the real and imaginary parts of the GOP, and, 

a„ = a„(Ap,Zp,Ar,E) , (20) 

fin = Pn{Ap,Zp,A„E) , (21) 

Jn = 0A5{^-^yAf+A]IUl). (22) 

Here, Ap, Zp, A,, and E are the mass number of the projectile, the proton number of the projectile, 
the mass number of the target, and the incident energy per nucleon, respectively. The "+1" in 
Eq. (I22l) is introduced to reflect the finite range of the A'^A'^ interaction. The factor (^) with 
n = 1 ~ 10 is taken to be an arithmetic progression from 0.5 to 1.0. 

We calculate the potentials for the incident even-even nuclei (^"^^C, ^^"^"^O, '^"''^Ne, ^"""^^Mg, 
22-48si^ 26-52s^ 30-62^^.^ ^^d ^^'^Ca) by the ^^C, 1^0, 28Si, 4°Ca 5^Ni, 90Zr, '^osn, and ^o^Pb targets at 
50, 60, 70, 80, 100, 120, 140, 160, 180, 200, 250, 300, 350, and 400 MeV/u. Here we only consider 
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a stable nucleus as a target. The desired optical potential is obtained by the spline interpolation 
over the incident energy and the target mass number. The parameters, a and /?, are available in 
Ref. [54]. The total number of the parameters provided is about 280,000. For a convenient use we 
also provide the program source "MGOP" in Ref. 
parameters and construct the optical potential. 



5411 which includes functions to generate the 
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FIG. 3: Comparison GOP with FMP and phenomenological optical potential (POP) for the ^^O + '"O 
system at £" = 70 MeV/u. The solid, dotted, and dashed curves are the FMP, POP, and GOP, respectively. 
The POP is taken from Ref je?] 



I6f 



Figure 12] shows the calculated FMP, phenomenological optical potential (POP), and GOP for 
the ^^O + '^O system at £ = 70 MeV/u. The GOP well fit to the FMP by a linear combination 
of 10-range Gaussian form, except for the most inner part of the real part {R = 0-2 fm). The real 
parts of FMP and GOP are slightly shallow in the comparison with that of POP. On the other hand, 
the FMP and GOP give the large imaginary potential. However, the both of the tail part almost 
becomes same strength. We will discuss this difference and the ambiguity of the imaginary part 
for the elastic cross section in the next section. 

In this paper, we do not provide the Coulomb part of GOP in a functional form and the standard 
Coulomb potential of a uniform charge is supposed to be used in the application of the nuclear part 
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1/3 1/3 

of GOPs proposed here. We use the radius of uniform charge, Rq = l-3(Ap + A/ ) fm. 



III. RESULTS 

A. Comparison with experimental data 

We now calculate the elastic scattering cross section with the folding model potential (FMP) 
and the present GOP fitted by 10-range Gaussian form. The results by these potentials are com- 
pared with the experimental data. The CEG07 folding model has only one parameter A''w Isoll and 
the optical potential is given by 



t/opt(i?) = Vf{R) + iN^^WpiR). 



(23) 



In this paper, we also discuss the energy and target dependences of this Nw value. 
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FIG. 4: Elastic scattering for the '^O + ^^C system at various energies. Abscissa q is the momentum transfer. 
The solid and dotted curves are the results by the FMP with the best fit and constant values (A'^i^ = 1.0), 
respectively. The dashed curves are the results by the present GOP with the Nw value. The experimental 



data is taken from Refs. 068 



TOD. 



We first test the FMP in the elastic scattering of some stable nuclei and investigate how the 
present FMP reproduces the experimental data free of renormalization factor. It should be empha- 
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size that we do not renormalize the real part of FMP in all the calculations shown below and we 
investigate the need of the renormalization to the imaginary part of the original FMP. 

Figures m |5] and |6] show elastic scattering cross section for the ^^O + ^^C, ^^O + '^O and '^C 
+ '^C systems at various energies. The solid and dotted curves are the results by the FMP with 
and without the renormalization factor for the imaginary part (Nw), respectively. The calculated 
cross sections reasonably reproduce the experimental data without the modification of Nw from 
unity, except for the low energy cases below 50 MeV/u. We thus conclude that the present FMP 
with CEGOVb gives a reasonable account of elastic scattering cross sections at the incident energy 
above 50 MeV/u in free of any adjustable parameters, although more perfect fit can be attained by 
a slight adjustment of Nw 

The dashed curves are the results by the present GOP with the same A^i^ values. The results 
by the present GOP at 20.6 MeV/u for 21.9 MeV/u for '^O+^^O and 25 MeV/u for '^C 

+ ^^C are not shown in these figures because we prepare the GOP in the energy range of 30 ~ 
400 MeV/u. The FMP and the present GOP give almost the same cross section, which implies 
that the fitting accuracy is enough to describe the elastic scattering given by the original FMP, 
although minor deviations from the original results are seen in high q region {q > 4 fm^) of the 
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FIG. 6: Same as Fig.|4]but for the '^C + ^^C system. The experimental data is taken from Refs. ll74l-l78h 



C + C scattering at higher incident energies above 100 MeV/u. The calculated elastic cross 
sections with GOP are close to the experimental data even without A^vi^ above 50 MeV/u and we 
decide to provide the GOP parameters for the incident energies above 50 MeV/u in the present 
paper. Here, we note that the FMP and GOP well reproduce the data despite the large difference 
of the imaginary potential seen in Fig. This implies that the experimental data cannot probe the 
imaginary-potential strength at short distance. 

Figure |7] shows the real part of the volume integral per nucleon pair of the GOP and POP for 
the '^O + '^C, '^O + '^O, and '^C + '^C systems. The volume integral per nucleon pair is defined 
as 

Jr*oo 
VGo?(R)R^dRIApAt. (24) 


The solid, dotted, and dashed curves are the results of GOPs for the '^O + '^C, '^O + '^O, and 
'^C + ^^C systems, respectively. The open circle, open square, and filled square are the results 
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FIG. 7: Volume integral per nucleon gair of die^GOP and POP for the '"O + '^C, ^"O + ^"O, and ^^C + ^^C 
systems. POPs are taken from Refs. 



74M78I1. 



of POPs for the ^^O + ^^C, ^''O + '^O, and ^^C + '^C systems, respectively. The ambiguity of 
the imaginary part is considered to be large as shown in Fig. |3] and mentioned in the previous 
paragraph. Therefore, the volume integral per nucleon pair only for the real part is shown in 
Fig. Ul The calculated volume integral per nucleon pair by GOP is consistent with the results 
by POPs. In addition, it is found that the volume integral per nucleon pair smoothly shifts from 
attractive to repulsive in the same way as shown in Ref. [|52|] . 

We also investigate the scattering by heavier target nuclei. Figures [8] and |9] show the elastic 
scattering of '^O and '^''Ar projectiles by various targets at E = 93.9 and 44 MeV/u, respectively. 
The FMP and the present GOP well reproduce the data up to backward angle without Nw, except 
for the '^O scattering by ^^Si case where a stronger absorption (Nw = 1 -6) is necessary to reproduce 
the data at backward angles. The origin of the exceptionally large value of Nw for the ^^Si target is 
unclear at present but the similar tendency was also reported in the case of a scattering by ^^Si [|49(| 
where the effect of strong collective excitation of the ^^Si target was suggested. 



B. Unstable nuclear region 



In this section, we compare with the elastic cross sections calculated by the Sao Paulo and 
Niigata densities. With the Sao Paulo and Niigata densities, the elastic scatterings of the Carbon 
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O elastic scattering 
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FIG. 8: Elastic scattering of the incident '^O particle for various targets at £" = 93.9 MeV/u. The solid and 
dotted cui^ves ai^e the results by the FMP with the best fit and constant values (Nw - 1-0), respectively. The 
dashed curves ai^e the results by the present GOP with the Nw value. The experimental data is taken from 
Ref. 



f. B. 



isotopes, Oxygen isotopes, and ^''Ne nucleus by the double-magic nuclear targets ('^O, '^''Ca, and 
2°^Pb) at 100, 200, 300, and 400 MeV/u are compared with the results by the present GOP. We fix 
the Nw value to be unity. 

Figure [10] shows the elastic scattering of the incident ^"^^^^C nuclei by the "^"Ca target at 100 
MeV/u. The dotted and solid curves are the results by the FMP with the Niigata and Sao Paulo 
densities, respectively. The dashed curves are the results by the present GOP. In this comparison, 
the FMP and GOP well reproduce the elastic scattering cross section calculated with the Niigata 
density, except for the result of the incident ^^C nucleus. The ^^C nucleus is considered and ob- 



served to have a two neutron halo structure [163 



8011 and the Sao Paulo density can not describe 



such exotic structure, that is why these densities give different results in the elastic cross sections. 
For the incident ^^C particle, the results of the GOP and FMP with Sao Paulo density overestimate 
that of the FMP with the Niigata density at backward angles. This reason comes from the differ- 
ence of the diffuseness of the densities. In this paper, we present the calculated results only in the 
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FIG. 9: Same as Fig. [8] but for the ^"Ar nucleus at E - 44 MeV/u. The experimental data is taken from 
Ref. B 



case by the ^"Ca target at 100 MeV/u. The others which are not shown are available in Ref. [|54ll . 

Figure [m shows the elastic scattering of the incident i^-^^Q nuclei by the '^"Ca target at 100 
MeV/u. For the Oxygen isotopes, the elastic cross section calculated with Sao Paulo density well 
reproduced the results with the Niigata one. In the case of the incident ^"^O particle, the difference 
is slightly seen in the cross section around backward angles. Because the neutron dripline nucleus 
of the Oxygen isotope (^"^O) is not very weakly binding system (B.E. ~ 4 MeV), that is, has no 
halo structure, the difference of the results between the Sao Paulo and Niigata densities is not as 
large as in the case of ^^C. 

Figure [I2] shows the elastic scattering of the incident ^"Ne particle by the "^^Ca target at various 
energies. Although the neutron number of the ^"Ne nucleus is a magic number, the shell gap 
vanishes around the ^"Ne nucleus and the region is known to be island of inversion [|8l|]. Therefore, 
the relation between the elastic cross sections with SP density and those with Niigata density 
of ^"Ne is similar to the ^"^O case. At 400 MeV/u, the result by Sao Paulo density reproduces 
that by Niigata density for the forward angles {q < 3), but the different results are given for the 
backward angles {q > 3). A similar result is obtained in the case of the incident ^"^O particle at 400 
MeV/u Jsj. 

Finally, Fig. [13] shows the elastic scattering of ™Ca by the "^"Ca target at various energies. 
According to Refs. [|82l- l85|] . we assume that the ™Ca nucleus is set to be a dripline nucleus. 
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FIG. 10: Comparison of the results by Niigata density with those by Sao Paulo (SP) density in the elastic 
scattering cross sections of the incident C isotopes by ^"Ca target at 100 MeV/u. 

Moreover, this nucleus is the heaviest projectile in this paper. For the heaviest projectile, the GOP 
well reproduce the elastic cross section calculated with FMP. Namely, it imply that the fitting 
accuracy of GOP is confirmed in the wide range region of the projectile. 

IV. CONCLUSION 

We have constructed a new global optical potential (GOP) in the framework of the CEG07 
folding model with Sao Paulo density. The incident particles are even-even isotopes of *^"^^C, 
i2-24q^ i6-38^g^ 22-4s^-^^ 26-523^ 30-62^^ 34-70(^a jj^g reliable energy range of 

GOP to be 50 - 400 MeV/u from the comparison with available experimental data. The present 
GOP can be used for heavy ion projectiles of stable and unstable nuclei over the wide range of 
nuclear chart with Z = 6 ~ 20, except for those on or near the driplines. The GOP parameters and 
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FIG. 11: Same as Fig.[TOl but for the incident O isotopes. 



the program source for constructing the GOP from the parameters are available in Ref. [|54ll . 
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